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lO ■ Abstract 

We prove a shape theorem for rotor-router aggregation on the comb, for a specific 
initial rotor configuration and clockwise rotor sequence for all vertices. Furthermore, 
as an application of rotor-router walks, we describe the harmonic measure of the 
rotor-router aggregate and related shapes, which is useful in the study of other 
growth models on the comb. We also identify the shape for which the harmonic 
measure is uniform. This gives the first known example where the rotor-router clus- 
ter has non-uniform harmonic measure, and grows with different speeds in different 
directions. 
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1 Introduction 

Rotor-router walks are deterministic analogues to random walks, which have been intro- 
duced into the physics literature under the name Eulerian walks by Priezzhev, D.Dhar 
et AL [PDDK96] as a model of self organized criticality, a concept established by Bak, 
Tang and Wiesenfeld [BTW88 . 
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In a rotor-router walk on a graph G, for each vertex x G G a cyclic ordering c(x) of 
its neighbours is chosen. At each vertex we have an arrow (rotor) pointing to one of 
the neighbours of the vertex. A particle performing a rotor-router walk carries out the 
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following procedure at each step: first it changes the rotor at its current position x to 
point to the next neighbour of x defined by the ordering c(x), and then the particle moves 
to the neighbour the rotor is now pointing at. 

The behaviour of rotor-router walks is in some respects remarkably close to that of random 
walks. See for example Cooper and Spencer [CS06] and Doerr and Friedrich 
[DF06] . 

In the present paper, we are interested in a process called rotor-router aggregation, defined 
as follows. Choose a root vertex o G G and let R\ = {o}. The sets R n are defined 
recursively, by 

R n+ i = R n U {z n } for n > 1, 

where z n is the first vertex outside R n that is visited by a rotor-router walk started at 
the origin o. The rotor configuration is not changed when a new particle is started at the 
origin. We will call the set R n the rotor-router cluster of n particles. 

Rotor-router aggregation on the Euclidean lattice Z d has been studied by Levine and 
Peres [LP09], who showed that the rotor-router cluster R n forms a ball in the usual 
Euclidean distance. On the homogeneous tree Landau and Levine |LL09j proved that, 
under certain conditions on the initial configuration of rotors, the rotor-router cluster R n 
forms a perfect ball with respect to the graph metric, whenever it has the right amount 
of particles. Kager and Levine |KL10] studied the shape of the rotor-router cluster on 
a modified two-dimensional lattice, which they call the layered square lattice. 

In each of these examples the fluctuations of the cluster around the limiting shape are 
much smaller in rotor-router aggregation than in the corresponding random growth model 
called internal diffusion limited aggregation (IDLA), where particles perform independent 
random walks before they settle and attach to the cluster. In the case of the homoge- 
neous tree and the layered square lattice, the fluctuations even vanish completely in the 
deterministic model. 

We will use the technique introduced in |KL10] in order to study rotor-router aggregation 
on the two-dimensional comb C2, which is the spanning tree of the two-dimensional lattice 
Z 2 , obtained by removing all horizontal edges of Z 2 except the ones on the x-axis. In other 
words, the graph C 2 can be constructed from a two-sided infinite path Z (the "backbone" 
of the comb), by attaching copies of Z (the "teeth") at every vertex of the backbone. 

We use the standard embedding of the comb into the two-dimensional Euclidean lattice 
Z 2 , and use Cartesian coordinates z = (x,y) G Z 2 to denote vertices of C2. The vertex 



(0, 0) will be the root vertex, see Figure 1(a) For functions g on the vertex set of C2 



we will often write g(x,y) instead of g(z), when z = (x,y). 

While C2 is a very simple graph, it has some remarkable properties. For example, the 
Einstein relation between the spectral-, walk- and fractal-dimension is violated on the 
comb, see Bertacchi |Ber06j . Peres and Krishnapur |KP04] showed that on C 2 
two independent simple random walks meet only finitely often; this is the so-called finite 
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■z = (x,y) 



(a) 



(b) 



Figure 1: (a) The two-dimensional comb C 2 - (b) The initial rotor configuration po 



collision property. 



The structure of this paper is as follows. In Section [2] we recall some preliminary results 
due to Kager and Levine [KL10] . which will be applied in order to prove the main 
result of the paper. In Section [3] we describe the shape of the rotor- router cluster on C 2 , 



for the initial rotor configuration po i n Figure 1 (b) Define 

B m = {(x,y) E C 2 : \x\ < m, \y\ < h(m — |x|)} for m G N, 



(1) 



for some function h : No — > No- The main result of the paper is the following. 

Theorem 1.1. Let R n be the rotor-router cluster of n particles on the comb C 2 , with 
initial rotor configuration as in Figure 1(b) and clockwise rotor sequence for all x € C 2 . 
Let B m be as in ([I]) with 

(x + l) 2 



h(x) 



Then, for all m > and n m = \B m \, the rotor-router cluster R Um satisfies R n 



B r , 



The main idea of the proof is to study rotor-router aggregation on the non-negative 
integers No, and then to glue different copies of No on the "backbone" oi C 2 , in order to get 
information on the behaviour of the rotor-router cluster. In the upcoming paper |HS11] 
the authors study IDLA on the comb. This gives another case where all fluctuations 
disappear in rotor-router aggregation when compared with IDLA. 

As an application of rotor-router walks, in Section H] we give a method to describe the 
harmonic measure of generic sets B m , of the form ([1]). For this, let 

dB m = {z£ B m : 3 neighbour y of z in C 2 , such that y ^ B m } (2) 

be the inner boundary of the set B m . The harmonic measure of B m is defined as the exit 
distribution of a simple random walk from the set B m . For z G dB m , we denote by v m ^{z) 
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the probability that a simple random walk started at the origin o G C2 hits dB m in z. In 
order to compute v mi0 (z), we consider a special rotor-router process, which allows us to 
obtain exact results in several cases. In particular, we identify the subsets of the comb 
for which the harmonic measure is uniform. 

Theorem 1.2. Let B m C C2 be as in ([T]) ; with h(x) = x 2 . Then the harmonic measure 
u m>0 of B m is the uniform measure on dB m . 

For the router-router cluster in Theorem 11.11 we are able to give asymptotics of the 
harmonic measure. We prove the following. 

Theorem 1.3. Let B m C C2 be as in ([I]) with h(x) = \_(x + 1) 2 /3J, and z = (z x ,z y ) G 
dB m . There exists a function e : No — > N with lim a; _ i , 00 = c, and < c < 1/2, such 
that for all m > the harmonic measure v m ,o{z) is proportional to e(m — \z x \). 

This gives the first example where the rotor-router cluster is not a set with uniform har- 
monic measure, and grows faster in the vertical direction than in the horizontal direction. 

We want to emphasize that it is not easy to apply this method in many cases, since it 
requires exact knowledge of the odometer function of the rotor-router walk, and at least 
some insight in the structure of the Abelian sandpile group of the set under consideration. 
The connection of the Abelian sandpile group to the rotor-router model has been estab- 
lished in the physics literature; see [PPS981 IPDDK96] . O ne can define a group based on 
the action of a particle which performs a rotor-router walk on the rotor configuration. 
This rotor-router group is abelian and isomorphic to the Abelian sandpile group. This 
isomorphism has been proven formally in |LL09] . For a self-contained introduction see 
the overview paper of Holroyd, Levine, et.al. [HLM+08] . 

2 Preliminaries 

Let (G, E{G)) be an infinite, undirected and connected graph, with vertex set G, equipped 
with a symmetric adjacency relation ~, which defines the set of edges E(G) (as a subset 
of G x G). We write (x,y) for the edge between the pair of neighbours x, y. In order to 
simplify the notation, instead of writing (G, E(G)) for a graph, we shall write only G, and 
it will be clear from the context whether we are considering edges or vertices. We denote 
by d(x) the degree of the vertex x, that is, the number of neighbours of x in G. Fix a 
nonempty subset S C G of vertices called sinks, and let G' = G\ S. 

The odometer function u(x) of the rotor-router aggregation is defined as the number of 
particles sent out by the vertex x during the creation of the rotor-router cluster R n of n 
particles. 

A rotor configuration on G is a function p : G' — > G, such that p(x) is a neighbour of 
x, for all x G G', that is (x, p{x)) G E{G). Hence, p assigns to every vertex one of its 
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neighbours. A rotor configuration p is called acyclic, if the subgraph of G spanned by the 
rotors contains no directed cycles. A particle configuration on G is a function o : G — > Z, 
with finite support. If o~{x) = m > 0, we say that there are m particles at vertex x. 
The rotor sequence at vertex x will be denoted by c(x) = (xo,xi, . . . ,Xd( x )-i) where all 
Xi ~ x and X{ ^ Xj for i ^ j, with i,j = 0, 1, ... , d(x) — 1. If y — Xi e c(x), for some 
i G {0, 1, . . . , - 1}, we denote by y + the vertex mod d(x) . 

Definition 2.1 (Toppling operator). Fix a vertex x G G' . For a rotor configuration p 
and a particle configuration a on G, we define the toppling operator F x , which sends one 
particle out of vertex x, by 

F x (p,a) = (p',a'), 
where the new rotor configuration p' is given by 

p\ y ) = [ p<yV ^ + i fy = x ' 

] p(y) otherwise, 
and the new particle configuration a' is given by 

(o-(y)-i tfy = x, 
a(y) + l ify = p'( x ), 
a(y) otherwise. 



So F x first changes the rotor configuration by rotating the arrow at x to its next position 
in the cyclic ordering c(x), and then it sends a particle to the vertex the rotor at x is now 
pointing at. The operation F x of toppling at some vertex x can be successful even if there 
is no particle at x. If this is the case, then a "virtual particle" is sent away from x and a 
"hole" is left there. If there is already a hole at x, the operator F x will increase its depth 
by one. In the normal rotor-router aggregation no holes are ever created during the whole 
process. A sequence of topplings {xk}k>i is called legal, if no holes are created when the 
vertices x^ are toppled in sequence. 

Note that the toppling operators commute, i.e., F x F y = F y F x for all x,y £ G'. This is 
the usual abelian property for rotor-router walks. While the final configuration is always 
the same, rearranging the order of the topplings can turn a legal toppling sequence into 
one that creates holes and virtual particles. 

Given a function u : G' — > N, let 

f u = n ^ 

x£G> 

where product means composition of the operators. Because of the abelian property, F u 
is well defined. 
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In order to prove a shape result for rotor-router aggregation on C 2 , f° r a specific initial 
configuration, we will apply a stronger version of the usual Abelian property of rotor- 
router walks, which has been recently introduced by Kager and Levine |KL10j . We 
state it here for completeness. 

Theorem 2.2 (Strong Abelian Property). Let po be a rotor configuration and o~o a particle 
configuration on G. Given two functions Ui, u 2 : G' — > N, write 

F Ui (p ,a ) = (puai), z = l,2. 

If o~i = a 2 on G' , and both pi and p 2 are acyclic, then u\ = u 2 . 

Note that the equality Ui = u 2 implies also that pi = p 2 , and moreover o\ = a 2 on all of G. 
This result allows us to drop the hypothesis of legality: each final particle configuration 
can only be achieved by an unique amount of topplings for each vertex, even if we allow 
virtual particles to be formed during the process. 

Friedrich and Levine [FLU] used the Strong Abelian Property to give an exact charac- 
terization of the odometer function of rotor-router aggregation. Recall that the odometer 
function u(x) at some vertex x represents the number of particles sent out by x during 
the creation of the rotor-router cluster. 

Theorem 2.3 (Friedrich, Levine). Let G be a finite or infinite directed graph, p an initial 
rotor configuration on G, and o"o = n ■ S a . Fix u± : G — > N, and let 

A* = {x e G : u*(x) > 0}. 

Further define p± and o~+ by 

Suppose the following properties hold 

(a) o±< I, 

(b) A* is finite, 

(c) <J*(x) — 1 for all x € A*, and 

(d) is acyclic on A*. 

Then is the rotor-router odometer function of n particles. 

Using Theorem l2.3[ for the proof of Theorem 1 1.1 1 it is enough to give an explicit formula for 
the corresponding odometer function, and to check if it satisfies all the required properties. 
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Figure 2: The first three fully symmetric configurations, consisting of n particles. The 
numbers on the arrows are the values of the odometer function u„. 



3 Rotor-Router Aggregation 



Consider now the rotor- router aggregation on the comb C2, and the initial rotor config- 



uration as in Figure 1(b) Through this section, B m will be the set defined in (CQ), with 
h(x) given by 

n(x) = 

Definition 3.1. Let (p, a) be the final configuration of the rotor-router aggregation process 
of \B m \ particles described in Theorem M . 1\ The configuration (p, a) is then called the m-th 
fully symmetric configuration. 



Figures [2] and |3] show examples of fully symmetric configurations. In Figure [31 one can 
observe that the fully symmetric configuration of \B 7 \ particles, as well as its corresponding 
rotor- router odometer function, are obtained by shifting "half" of the configuration of \B§\ 
particles one step in the direction of the positive resp. negative x-axis and filling in the 
values for the "tooth" corresponding to x — 0. It turns out that this is true for all fully 
symmetric configurations (with the exception of the first 3). This property will play an 
important role in the proof of Theorem 11.11 

For proving Theorem 11.11 an exact expression for the cardinality of the sets B m is needed. 



Lemma 3.2. Let B m be the set defined in ([T]) ; with h(x) 
the cardinality of B m is given by 



Qr+l) 2 



. Then, for allm > 



\B m \ = - [4m 3 + 12m 2 + 24m + 5 + 2((m + 2) mod 3)] . 
9 



(3) 



Proof. In order to simplify the statement of the Lemma, we have to distinguish three 



7 



cases, namely for m = 3k + i, with i = 0,1, 2. The right-hand side of (j3J) is then equal to 
N (k) = 12k 3 + 12k 2 + 8k + l, for m = 3k 

iVi(Jfe) = 12k 3 + 24k 2 + 20k + 5, for m = 3k + 1 (4) 
N 2 (k) = 12k 3 + 36k 2 + 40k + 15, for m = 3k + 2. 

Moreover, in the three cases m = 3k + i, with i = 0,1, 2, the function h can be written as 
follows. 

h(3k) = 3k 2 + 2k 
h{3k + 1) = 3k 2 + 4k + 1 (5) 
/i(3A; + 2) = 3/fc 2 + 6A; + 3. 

We prove by induction on fc. The base case m = 1 is immediate from the definition of 
B m . With 

|Snn-il = |-B m |+2[/i(m)+/i(m + l) + l] 
follows the inductive step. □ 

Next, we will find an exact formula for the odometer function of the rotor- router aggre- 
gation defined in Theorem and we shall prove its correctness using Theorem 12.31 For 
this, we first have a detailed look at the rotor-router process on the non-negative integers. 



3.1 Rotor-Router on the non-negative Integers 

For a better understanding of the rotor- router process on the comb C 2 , we first analyse it 
on the half- line, where it is very simple. Consider G = Nq, with sink vertex 0, and the 
initial rotor configuration p : N — > Nq given by 

po(y) = y + 1) for a11 y e n. 

Let Ri = {1}, and define a modified rotor-router aggregation process R n recursively as 
follows. Start a rotor-router walk in 1, and stop the particle when it either reaches the 
sink 0, or exits the previous cluster R n -\. Denote by z n the vertex where the n-th particle 
stops, and by p n and u n the rotor configuration and odometer function at that time. 
Then, 

R n -i U {z n }, if Zn ^ 
Rn-i, otherwise. 



Rn 



Obviously R n = {1, . . . , h(n)} for some sequence h(n). Since po is acyclic, all rotor 
configurations p n are acyclic and have the form 



Pn(y) 




< y < f(n) 
otherwise, 
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n = 161 



n = 237 



Figure 3: The 6 th and 7 th fully symmetric configurations, consisting of \B$\ and \B 7 \ 
particles. The numbers are the values of the odometer function u 6 and u 7 , respectively. 
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Figure 4: The first steps of the process R n on N. The dots mark the vertex where the 
current particle stopped. 



for some numbers < r(n) < h(n). Here, r(n) represents the vertex where the rotors 
change direction: all rotors up to r(n) point inwards (J,), and all rotors from f(n) + 1 up 
to h(n) point outwards (f). 

For numbers h, r and y in No, with < r < h define the function u as 

{f(h-y) + e(r-y), l<y<r 
f(h-y), r<y<h (6) 

0, otherwise, 

where the functions e and / are given by e(y) — 2y + 1 and f(y) = y(y + 1). 

The odometer function u{ri) of the rotor-router process R n can now be defined in terms 
of (JBj) by setting 

h(n) = max jfc 6 N : ^+ ^ < n | , f(n) = n - (7 ) 

and 

"n(y) = u(h(n),f(n),y), (8) 

for /i(n) and f(n) as defined in (j7]). It is easy to verify by induction that u n correctly 
describes the odometer function of R n . See Figure H] for a graphical representation of the 
process R n . 



3.2 Rotor-Router on the Comb 

Since the rotor-router aggregation on the "half-teeth" of C2 behaves like the process R n 
from the previous section, it is enough to determine the numbers h and r in (J6]), which now 
depend on x and the number of particles, in order to fully specify the odometer function 
on C 2 for points off the x-axis. 
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Let B m as defined in (P), with h(x) 



and define r(x) by 



x£{0,l} 

r(x) = <j ±(x 2 - 7x + 10), i = 2 mod 3 



±(x 2 - x + 6), 



otherwise. 



Define u m : £? m — > N by 



where 



u'(x,y) 



u m (x,y) = u'(m - \x\, \y\), 

u(h(x),r(x),y), 
2f(h(x))+2e(r(x))-2-t {x=2} , y = 0, 



y>0 



with m as in ([6]), and r(x) defined as above, and e(x) = 2x + 1 and f(x) = 
We claim that u m is the odometer function for rotor-router aggregation of \B n 
on the comb. 



(9) 



(10) 



(11) 

x(x + 1). 
j particles 



3.2.1 Proof of Theorem [TTT1 



Let _B m as defined in (CQ), with h(x) = [(x + l) 2 / 3J as in Theorem ll.il and u m defined as 
in (HUD. 



Lemma 3.3. Lei po ^ e the initial rotor configuration defined in Figure 1(b), and the initial 
particle configuration o"o = \B m \ ■ 5 . Furthermore, define p m and o~ m as 

(p m ,o- m ) = F Um (p ,a ), 

with u m as in (|T0|) . If a clockwise rotor sequence is assumed for all vertices, then u m is 
the odometer function of the rotor-router aggregation with \B m \ particles and moreover 
cr m = lfl m , for all m > 3. 



Proof. To verify that u m is indeed the odometer function of this rotor-router process, we 
need to check the four properties of Theorem 12.31 with 

A* = B m \ OB 

Til ■ 

where dB m is the inner boundary of B m defined in ([2]). The set A* is obviously finite. 
For those vertices z £ A* that have neighbours in B m \ A*, we have by (jHJ) and ([6]) that 
u m {z) < 3 if z is not on the x-axis, and u m (z) = 4 otherwise. In both cases at most one 
particle is sent to some vertex outside of A+, hence <J m {z) < 1 for all z £" A+. 

Next we verify that the final particle configuration a m is equal to 1 on A*. Since by 
definition u m is symmetric, it is enough to consider only one quadrant. Additionally, we 
shift the coordinate system such that the point (— m, 0) lies at the origin, which means 
that we can work with the function u', defined in ( fTTT) . Since u' does not depend on the 
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parameter m, most of what follows holds independently of m. Only for the center (m, 0) 
of the set B m (Case 4), we need to take the parameter m into account. Let z = (x, y) G A* 
with x,y > 0. We distinguish several cases. 

Case 1. y > 2: For vertices (x,y), with y > 2, the rotor-router aggregation behaves 
as the process defined in Section 13.11 From Figure HI and due to the fact that the 
final rotor configuration restricted to each "tooth" is acyclic, there are again four possible 
situations: 

(a) The rotors at the vertices (x, y — 1), (x, y) and (x, y + 1) all point outwards (t). This 
is the case when r(x) < y — 1, hence the vertex (x, y) receives \u m {y — 1) + \u m (y + 1) 
particles from its upper and lower neighbours, and it sends u m {y) particles. That is, 

<r m (x, V) = \ [f(h(x) -y + l) + f(h(x) - y - l)] - f(h(x) -y) = l. 

(b) The rotors at the vertices (x, y — 1), (x, y) and (x, y + 1) all point inwards (|). Hence 
r(x) > y + 1 and, comparing the numbers of incoming and outgoing particles, we have 

a m (x,y) = ~[f(h(x) -y+l)+ e(r(x) - y + l) + f(h(x) -y-l) 
+ e(r(x) - y- l)] - f(h(x) - y) - e(r{x) - y) = 1. 

(c) When the rotors from 1 to y — 1 point inwards (J,) and from y to h(x) point outwards 
(t), then r(x) — y — 1, and we have 

(?m(x,y) = ^[f{h{ x ), r ( x ),y- 1) + e (°) - i] 

+ ^/( /i ( s ) ) r (a ; ) ) l/ + l) - f{h(x),r(x),y) = 1. 

(d) The last case which can appear is when all rotors from y to point inwards (J,), and 
from y + 1 to h(x) outwards (f). Then r(x) = y and 

X 

a m (x,y) = ~[f(h(x),r(x),y-l) + e(l) - l] 

+ ^/( /i (^)> r (a ; ) ) z/ + 1) - f{H x ), r { x )>y) -1 = 1. 

Therefore a(x,y) = 1 for all (x, i/) G -B m with y>2. Moreover, no closed cycle is formed 
by these vertices in the final rotor configuration p m . 

Case 2. y = 1: Consider a m (z) for the vertex z = {x, 1). For x > 9 the number of 
inwards pointing rotors r(x) is always greater than 2. So with the exception of a finite 
number of exceptional points (x G {1, 2, 5, 8}), all relevant rotors on the teeth are pointing 
inwards (l) and the vertex z receives [|m'(x,2)] particles from its upper neighbour. For 
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x > 3, the number u'(x, 0) is divisible by 4, so all neighbours of (x, 0) receive exactly the 
same amount of particles. Hence 

a m (x, 1) = ±u'(x, 0) + ±(u'(x, 2) + 1) - u'(x, 1) = 1 [2/(/i(x)) + 2e(r(x)) - 2] 

+ ^ [/(%) - 2) + e(r(x) - 2) + 1] - f(h(x) - l) - e(r(x) - l) 

= ~ [2/i(x)(/i(x) + 1) + 2(2r(x) + 1) - 2] 

+ ^[(h(x) - 2){h{x) - 1) + 2{r{x) - 2) + 2] 
- (%) - - 2(r(x) - 1) - 1 = 1. 

if z is non-exceptional. At the exceptional points z = (x, 1), for x e {1,2,5,8}, the 
correctness of the function w' can be verified by direct computation. 

Case 3. x ^ m and y = 0: On the x-axis, the points z = (x, 0) for x G {2, 5} are 
again exceptional and need to be checked separately. The case x = does not have to be 
checked at all, and x = 1 has already been checked at the start of the proof. 

For x {2, 5}, the vertex z = (x, 0) receives particles from (x — 1,0), (x + 1,0), (x, 1), 
(x, —1). Here u'(x — 1,0) and u'{x + 1,0) are again both divisible by 4. By symmetry 
u'(x, 1) = u'(x, —1), and the number of inward pointing arrows r(x) > 1 in this case, hence 
z receives u'(x, 1) + 1 particles from its upper and lower neighbours combined. Thus 

cr m (x, 0) = -u'(x - 1, 0) + ^u'(x + 1, 0) + u'(x, 1) + 1 - u'(x, 0) 

= -[2f(h(x-l)) +2e(r(x- 1)) -2} + - [2f(h(x + 1)) + 2e(r(x + 1)) - 2] 
+ f(h(x) - 1) +e(r(x) - 1)) + 1 - [2f (h(x)) + 2e(r(x)) -2]. 
Using that f(x) = x(x + 1) and e(x) = 2x + 1 we get 

a m (x, 0) = - [h 2 {x - 1) + h 2 {x + 1) - 2h 2 (x)] + - [h{x - 1) + h{x + 1) - 6h(x)] 

2 2 (12) 

+ [r(x — 1) + r(x + 1) — 2r(x)] . 
In order to check a m (x,0) = 1, we have to substitute in equation ( [121 the function 



h(x) 



(x+l) 2 



and the corresponding branch of the function r(x) given in equation (J9j) 



3 

depending on the congruence class mod 3 of x. We have to check all three cases separately. 
In all cases <J m (x, 0) = 1 holds. 

Case 4. Midpoint z = (m,0): Everything until now was independent of the number 
of particles \B m \. Since u m is created from u' by translation and reflection, the vertex 
z = (m, 0) after translation corresponds to the origin of the cluster. At the beginning of 
the process, \B m \ particles are present at z, so (Jq(z) = \B m \. We assume that m is big 
enough, so that none of the neighbours of z is an exceptional point. 
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By symmetry, z = (m, 0) receives ~u'(m — 1,0) particles from its neighbours on the x 
axis, and u'(m, 1) + 1 particles from its neighbours on the teeth. Hence 

cr m (m, 0) = cr (m, 0) + -u'{m — 1,0) + u'(m, 1) + 1 — u'(m, 0) 

= \B m \ + l - [2f(h(m - 1)) + 2e(r(m - 1)) - 2] + f(h(m) - l) 

+ e(r(m) - l) - 2f(h(m)) - 2e(r(m)) + 3. 

Here one has to check again each congruence class mod 3 separately Substituting the 
formulas for \B m \ obtained in (j3]), into the previous equation, gives the desired result 
a m (z) = 1. 

Finally, we need to check that the final rotor configuration p m is acyclic. We work again 
with shifted coordinates. It is clear from the previous section that p m restricted to each 
"tooth" is acyclic. Hence it suffices to check that no cycles are created by rotors on the 
x-axis. If z = (x,0), the odometer u m (z) is divisible by 4, except when x = 2. So the 
rotors at these vertices point in the same direction as in the initial configuration p Q . The 
odometer at the exceptional point w = (2, 0) is u'(w) = 23 = 3 (mod 4) independent 
of m. Hence, this rotor points in the direction of one "tooth". If the rotor at position 
(2,1) points towards the x-axis, it creates a directed cycle. By fl9]), we have r(2) = 0, 
which means that all arrows on this "tooth" are pointing outwards. Hence the rotor at w 
does not close a cycle. See Figure [3] for a visualisation of the rotor configurations under 
consideration. 

Therefore all properties of Theorem 12.31 are satisfied and this proves the statement. □ 

Proof of Theorem \l.l\ In the case m < 2, the statement of the Theorem follows by direct 
calculation of the respective aggregation clusters, see Figure EJ For m > 3 it follows from 
the previous Lemma. □ 



4 Harmonic Measure 

In this section, as a direct application of rotor-router walks, we compute the harmonic 
measure of the generic set B m C C2 defined in ([[]). The harmonic measure of B m is the 
hitting distribution of the set dB m for a simple random walk on C2 starting at the origin 
o. 

We shall first describe the method for finite subsets B of general graphs G, and then we 
apply it to the case of the comb C 2 and subsets B m of the type defined in (Q. In Theorem 
ll.2[ we identify the shape for which the harmonic measure is uniform. We point out 
that this shape does not coincide with the rotor-router aggregation cluster from Theorem 
11.11 We will also describe the asymptotics of the harmonic measure for the rotor-router 
aggregation clusters. 
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In order to estimate the harmonic measure, we shall use an idea of Holroyd and Propp 
[HP 10] . which they used to show a variety of inequalities concerning rotor-router walks 
and random walks. The method assigns a weight to the particle and rotor configuration 
of a rotor-router process, which is invariant under routing of particles in the system. 



4.1 Rotor Weights 

Let G be a locally finite and connected graph. Start with a particle configuration er : 
G — > Z and a rotor configuration p : G — > G such that po(x) = xq for all x e G, that is, 
all initial rotors point to the first neighbour in the rotor sequence c(x). We further assume 
that o"o has finite support, i.e., there are only finitely many particles in the system, so 
that we don't need to deal with questions of convergence. We will route particles in the 
system, and this gives rise to a sequence (p t , o~t)t>o of particle and rotor configurations at 
every time t. To each of the possible states (p t , a t ) of the system, we will assign a weight. 

Fix a function xjj : G — > R. We define the particle weights at time t to be 

W P (t) = 5^ < 7 t (x)V(x). (13) 

xeG 

Further define the rotor weights of vertices x G G as 

, n f0, forA; = 

w(x,k) = < * . (14) 

[w(x, k - 1) + ip(x) - ip{x kmod d(x)), for k> 0, 

where Xi is the i-th neighbour of x in the rotor sequence c(x). Notice that, for k > d(x), 

w(x, k) — w(x,k — d(x)) — d(x) A ip(x). (15) 

Here Aij}(x) represents the Laplace operator which is defined as 

The total rotor weights at time t are given by 

W R (t) = ^u;(x,^(x)), 

where u t (x) is the odometer function of this process, that is, the number of particles sent 
out by the vertex x in the first t steps. Note that po is chosen in such a way that for all 
t > and x € G, if i = Mt(x) mod g?(;e), then Xj = pt(x). 

It is easy to check that the sum of particle and rotor weights are invariant under routing 
of particles, i.e., for all times t, t' > 

W P (i) + W R (t) = W P (f) + W R (f). (16) 
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4.2 Harmonic Measure for finite subsets of graphs 

As before, let G be a locally finite, connected graph, and let B be some finite subset of 
G. Write 

dB = [x G B : 3y G' B with x ~ y} 

for the inner boundary of B, and B° = B\ dB. The vertices of (95 will represent the sink 
S. 

Similarly to Definition 12.11 of Section [2j we define the particle addition operator E x , for 
each vertex x G B°, as follows: for a rotor configuration p, let 

EM = p\ 

where p' is the rotor configuration obtained from p by adding a new particle at vertex x, 
and letting it perform a rotor-router walk until the particle reaches a vertex in dB for 
the first time. By the abelian property of rotor-router walks the operators E x commute, 
and they can be used to define an abelian group, see |HLM+08] for details and [HLM+081 
Lemma 3.10] for the proof of the following statement. 

Lemma 4.1. The particle addition operator E x is a permutation on the set of acyclic 
rotor configurations on B° . 

The rotor-router group of B° is defined as the subgroup of permutations of oriented 
spanning trees rooted at the sink (that is, acyclic rotor configurations) generated by 
{E x : x E B°}. For every finite graph B° the rotor-router group is a finite abelian group, 
which is isomorphic to the abelian sandpile group. See once again [HLM+08] for details. 

Consider the simple random walk (X t ) t > on G, i.e., a Markov chain with state space G, 
and transition probabilities given by 

p(x, y) = , for all x,y E G, with x ~ y. 
d(x) 

Then X t is a G- valued random variable, and represents the random position of the random 
walker at the discrete time t. For given x G G, we write for the law of a random walk 
starting at x. Consider the stopping time 

T = inf {t > : X t G dB}. 

For zedB, let 

vM = ^{Xt = z], 

be the harmonic measure at z with starting point x, that is, the probability that a random 
walk starting at x hits dB for the first time in z. 

Take the harmonic measure itself as the weight function. More explicitly, fix a vertex 
z G dB, and define the weight function ip(x) as 

ip(x) = ip z (x) = v x {z). 
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Let us define the following process. Start with n particles at the origin o, and an arbitrary 
acyclic rotor configuration p . Let the particles perform rotor- router walks until they reach 
a vertex in dB for the first time, where they stop. Denote by t* = t*(n) the number of 
steps this process takes to complete, and for each w G dB, write e(w) for the number of 
particles that stopped in w at the end of this procedure. We denote by u the normalized 
rotor- router odometer function of this process, that is, for all x G B, 

number of particles sent out by x 



Using the invariance of the sum of rotor and particle weights under rotor-router walks, as 
in (TIB]) , we get 

nif>(o) = < w )^( w ) + W R (t*), (17) 

wEdB 

since W R (0) = 0, W P (0) = nif>(Q) and W P (t*) = J2 w&9 b e{w)i/>(w). Equation 07jl 
reduces to 

niP(o) = e(*)+W R (t*), (18) 
because ip(w) = v w {z) = 6 w (z), if w G dB. 

The initial rotor configuration p is chosen to be acyclic. Therefore, there exists a number 
n such that, after all n particles performed their rotor- router walks, all rotors in B° made 
only full turns, i.e., p$ = p*. This claim follows from Lemma [4.11 Hence, n is a multiple 
of the order of E Q in the rotor-router group. Since ip is a harmonic function on B°, using 
a n with the above property gives Wa(f) = 0, which together with ( 118)) leads to 

n-u (z) = e(z). (19) 

Thus, the harmonic measure v Q of B is proportional to the number of particles which 
stopped at dB. While a number n with the right property is difficult to calculate, we can 
still use equation (fl9~|) in order to derive asymptotics of the harmonic measure of subsets 
B m of the comb C2, and in some cases even to calculate it explicitly. 



4.3 Subsets of the Comb 

Let us consider subsets B m of C 2 > of the type defined in (JT]), with generic positive function 
h : No — > No- Recall here the definition of B m . 

B m = { (x, y) G C2 : |x| < m, \y\ < h(m — \x\)} for m G N. 

By construction, all rotors make only full turns if we perform the rotor-router process from 
Section \A.2\ for the set B m . This implies that the corresponding normalized odometer 
function u is harmonic outside the origin and its Laplacian is given by 

A4(„) = <- t " eB ™\( aB ™ U {°» (20) 
w — o, 
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and u(w) = 0, for w G dB m . Write u m>0 (w) for the harmonic measure of B m and e m {w) 
for the number of rotor- router particles stopped in w G dB m . By symmetry of the set 
B m , it is clear that also e m (w) and z/ m>0 (u>) are symmetric. More precisely, if w = (x,y) 
and w' = (\x\, \y\) then 

e m (w) = e m (w') and v m , {w) = v m ,o{ w ')- 

Hence it is enough to work in one quadrant. We will choose the second quadrant (i.e. 
x < and y > 0), and for simplicity of notation shift the set B m by m in the direction of 
the positive x axis, such that the leftmost point of B m has coordinate (0, 0) and its center 
o has coordinate (m, 0). So, the set under consideration is now 

B° = {(x,y) GC 2 : < x < m, 0<y< h(x)} 

We will also use u for the normalized odometer function and e m for the number of particles 
which hit boundary points in the shifted coordinate system. Additionally, since e m is 
defined only on dB m we write e m (x) = e m (x, h(x)), for < x < m. 

Solving the Dirichlet problem ([20]) on the "teeth" of the comb, gives for (x, y) G B^, 

u(x,y) = e m (x) ■ (h(x) -y). (21) 

On the x axis, for (x, 0) 7^ o, the harmonicity gives 

u(x + 1,0)+ u(x -1,0) + 2u(x, 1) = Au{x, 0), 

which together with (1211) leads to the following recursion for e m (x) and < x < m: 

e m (x + l)h(x + 1) + e m (x - l)h(x - 1) - 2e m (x) (h{x) + l) = 0. (22) 

We are now ready prove Theorem II. 2\ which characterizes sets of uniform harmonic mea- 
sure. 

Proof of Theorem M.B. Let B m be defined as in ([1]), with h(x) = x 2 . From fT22|) we get the 
recursion 

e m (x + l)(x + l) 2 + e m (x-l)(x-l) 2 -2e m (a;)(x 2 + l) = 0, for < x < m. (23) 

Since h(l) = 1, the vertex z — (1, 0) in the shifted coordinate system has three neighbours 
on the boundary dB m . By construction, the rotor at (1, 0) makes a number of full turns, 
hence all of these three neighbours receive the same amount of particles from z. Therefore 
e m (0) = e m (l). By induction, it is easy to see that the sequence e m (x) is constant. 
Assuming e m (x — 1) = e m (x), the recursion fl23|) reduces to 

e m (x + l)(x + l) 2 - e m (x)(x + l) 2 = 0. 

which implies that e m (x + 1) = e m (x). Because e m (x) is by construction proportional to 
the harmonic measure z/ m o , we get the claim. □ 

In general we can compute the harmonic measure for all sets B m , where ( !24l) can be solved 
exclicitly. 
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4.4 Harmonic Measure of the Rotor-Router Cluster 



This section is dedicated to the proof of Theorem II .31 Recall that the rotor-router cluster 
obtained in Theorem 11.11 is a set of type B m , as defined in ([I]), with h(x) = \_(x + l) 2 /3j. 

Like before, let e m (x) be the number of particles stopped at boundary points (x, h(x)) G 
dB m in the rotor-router process defined in Section I4.2L By linearity also the normalized 
sequence e(x) = f^^y is a solution of the recurrence ff22l) . and since h(l) = 1 we have 
e(0) = e(l) = 1. Hence the function e(x) is independent of m. Rewriting (|22l) in this 
case, we get a linear recurrence with non-polynomial coefficients. While an explicit answer 
is not feasible, we can derive asymptotics of the special solution e(x), by converting 
the recurrence into an equivalent system of linear differential equations. We will prove 
Theorem 11.31 by showing that the function e(x) has linear growth, which is accomplished 
in the next two Lemmas. 

Lemma 4.2. There exists a constant c < \ such that 

e(x) 

> c as x — >• oo. 

x 

Proof. Substitute e(x) = for x > 0, which transforms (122]) into 

e(x - l)(x - l)h(x - 1) + e{x + l)(x + l)h(x + 1) - 2e{x)x(h{x) + l) = 0. (24) 

The sequence e(x) converges if and only if e(x) grows at most linearly. Since e(x) is 
positive by construction, it suffices to check that e(x) is decreasing. For this, consider the 
auxiliary function h'(x) = -^4p |. We have to distinguish three cases 

{h'(x), x = mod 3 

h'(x), x = 1 mod 3 

h'(x) + |, x = 2 mod 3 

We prove the monotonicity of e(x) by induction. Assuming e(x) < e(x — 1) for x = 
mod 3 we show that e(x + 3) < e(x + 2) < e(x + 1) < e(x). The induction base follows 
by calculating the first elements of the sequence. 

Case 1. Assume x = mod 3 and e(x) < e(x — 1). Then (|24j) can be rewritten as 
e(x + l)(x + l)h'(x + 1) = 2e(x)x(h'(x) + l) - e{x — l)(ar — 1) (h'{x - 1) + |) . 
Using the induction hypothesis and the definition of h'(x), we get 

e(x + l)<f (x)-e(x), (25) 
with fo(x) = x 2^x+i < w hich implies e(x + 1) < e(x). 
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Case 2. Assume x = 1 mod 3 and e(x) < fo(x — 1) • e(x — 1). Like before, rewrite (T2 
as 

e(x + l){x + l)(h'(x + 1) + |) = 2e(z)x(/i'(x) + l) - e(x - l)(x - l)h\x - 1). 
This gives, by ff25l) 

e(x+l)</ 1 (a;)-e(a;) J (26) 

for 

~ 2x(h'(x) + 1) - /o(x - l)- 1 ^ - l)h'(x - 1) 

H{X) (x + l)(h'(x + l) + l) 

x 2 + 3x 
~ x 2 + 3x + 2 < ' 

which implies e(x + 1) < e(x). 

Case 3. Finally, assuming x = 2 mod 3 and e(x) < fi(x — 1) • e(x — 1), we get 
e(ar + l)(x + + 1) = 2e{x)x(h'{x) + |) - e(x - l)(z - - 1). 

Applying (126]) . we obtain 



e(x + l)<f 2 (x)-e(x), (27) 
for the function f 2 (x) = J^X^X^Ze < 

This shows that e(a?) is decreasing and therefore convergent, which also means that there 
exists a constant c such that — > c, as x — > oo. The fact that c < \ follows by 
computing the first few values of the sequence e(x), using e(0) = e(l) = 1 as starting 
values in the recursion (1221) . By monotonicity we then get e(x) < | for all x > 20. □ 

The next result shows that the function e(x) has at least linear growth. 
Lemma 4.3. The constant c in Lemma \4-S\ is strictly positive. 

Proof. To show that c > 0, we use singularity analysis of linear differential equations. For 
this, we split e(x) into three sequences modulo 3, i.e., for k G N write 

Gi (k) = e{3k + i) for i = 0,1, 2, 

and rewrite ( 122]) for each congruence class of x mod 3 in terms of k. This leads to a system 
of linear recursions which can be written in matrix form as 

A k -e(k-l) = B k -e(k), (28) 
with e(k) = (e (k), ei(k), e 2 (fc))', and the matrices A k and B k given as 

/0 3fc 2 - 2k -6k 2 - 2\ / -3k 2 -2k \ 



3fc 2 

Vo y 



Blr. 



6fe 2 + 4fc + 2 -3fc 2 -4fc-l 
V 3fc 2 + 2k -6k 2 -8k -4 3k 2 + 6k + 3/ 
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The initial values are given by e(0) = (1,1, |)*. Denote by E{(z) = ^2 k>0 ei(k)z k the 
generating function of ei(k), i = 0, 1, 2. Using the identities 

d 

~^^kei(k)z k = z—Ei(z) and k 2 ej(k)z 

k>0 ' k>0 



z 



the matrix recursion ( 12811 can be transformed into the following system of linear differential 
equations for the generating functions Ei(k) 

C ■ E{z) = b, (29) 

where E(z) = (E (z) , Bk(z) , E 2 (z)) i , and C is a matrix of linear differential operators 
given as 



C- 



5 9 3z — 
dz dz 2 



102- 



_3 

'al 



V 



9 , a 2 

5z—+3z 2 — 
oz oz z 



a , a 2 
a , a 2 

1 + 72— +32 2 -^ 

az oz* 

a , a 2 

_4_14 Z _ 6z 2 



182- 



9 9 2 \ 
62 2 -> 

82 a2 2 
, a , a 2 

32 + 92 2 — + 32 3 — 

02 02 z 

a , a 2 
3 + 9z aI + 32 2 a? / 



ei(0) -2e (0) 




32 dz 2 

To solve ( 129]) asymptotically, we consider C as a matrix with entries in the Weyl algebra, 
that is, the noncommutative ring of linear differential operators with polynomial coeffi- 
cients, see |Lam91] . We can perform a division-free Gauss elimination over this ring to 
transform C into row echelon form, which gives a single differential equation only involv- 
ing E2(z). The actual computations were performed using the computer algebra system 
iCAf]. 
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The result is a differential equation of order 7 for E 2 (z): 



f(z + 2)(z 



+!?>(* -1)V (z 2 + z-§)^ 



127725 



(z-lf 



+31785(z- l)z 2 (; 



I 27531 
8 

Z* ■ (z 4 

_ 100697 
42380 



(z- 1)V (r 5 



24 2 _ 7149 , 3826\ cfi_ 
437 3059 ~ t ~ 3059/ dz 5 



z 4 - 



50039 3 
42575 

1164 



82401 r 2 I 132307 _ 



42575 



42575 



3 I 36215 2 _ 5651 

10595 ' 8476 1630 



+239702 • (i 



117579 ^5 _|_ 



114057 ^4 _|_ 15053 Ji 



31960 * 1 31960 ~ 6392 ~ 31960 

3 



208329^2 _|_ 59229 ^ 



15980' 



38554 \ 8 4 
42575/ dz 4 

I 6234 \ 8 3 
10595/ 32 s 

1243 \ a 2 

3995/ 32 2 



/LOQP; (+6 _ 1354 5 I 1843 r 4 _ 4479 

-±t)00 329 4, t * 3290 



12209 2 , 1466 
3290 329 



32 \ _a_ 

329/ 32 



+ 105 (i 



494 4 _|_ 881„3 _ 201 2 , 4411 
105 ' 105 70 ~ r 210 



1006 N 
105 / 







(30) 



Using singularity analysis for linear differential equations, we can derive asympotics of 
e 2 (k). See Flajolet and Sedgewick |FS09] Theorem VII. 10] for details. The coeffi- 
cient of the highest order term is given by 

^-{z + 2)(z-l) 5 z 6 , 



^ttp : / /f ricas . sourcef orge . net 
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hence the dominant non-zero singularity £ is equal to 1. Since all coefficients in ( 130]) are 
given in factorized form, it is immediate that £ is a regular singularity. Calculating the 
indicia! 'polynomial for the singularity £ gives 

10) = e 7 - 176 6 + 996 5 - me 4 - 22oe 3 + wue 2 - 7209. 

For the definition of a regular singularity and the indicial polynomial, see once again 
Flajolet and Sedgewick [FS091 Chapter VII.9]. The roots of 1^9) are -2, 0, 1, 3, 4, 5 
and 6. Since they differ by integers, the asymptotics of e2(k) is given by 

e 2 (k) ~ c-C k k l3 log l k, 

where I is an integer and (3 is the biggest solution of the equation J(— /3 — 1) = 0, see 
[FS09j page 521, equation 118]. The ~ sign means "approximately equal" (in the precise 
sense that the ratio of both terms tends to 1 as k gets large). 

In our case (3 = 1, and we have 

e 2 (Jfe) ~ c- k\og l k, (31) 
and this proves the desired. □ 

While it is not known how to calculate the constant I in (1311) in the general case, from 
Lemma [4.21 we already know that C2{k) grows at most linearly, hence / = 0. Therefore, 
Lemma 14.21 and Lemma 14.31 together imply Theorem 11.31 
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